Decoherence-assisted transport and quantum criticalities 
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We study the dynamics of a two-level quantum system interacting with an external environment 
that takes the form of an XY spin chain in the presence of an external magnetic field. While 
the presence of the bath itself can enhance the transition probability from the lower level to the 
upper level of the system, we show that this noise-assisted phenomenon is sensitive to a change 
of the quantum phase of the environment. The derivative of the transition probability displays a 
maximum in correspondence with the critical value of the applied field both in the case of isotropic 
and anisotropic chains. 
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The unavoidable interaction of any realistic quantum 
system with its environment is responsible for decoher- 
ence and dissipation, which are usually considered detri- 
mental effects for the implementation of quantum infor- 
mation and quantum communication tasks 0, Q ■ How- 
ever, recent experimental studies have shown that some 
photosynthetic complexes exhibit long excitonic coher- 
ence despite the interaction with their own environ- 
ment 3 5]. Furthermore, the possibility of decoherence- 
assisted quantum transport was recently discovered in a 
series of theoretical Hamiltonian models. Therefore, fully 
understanding the role played by the noise in assisting 
quantum transport phenomena has assumed a crucial im- 
portance and has become a rich subject of investigation 
(see, for instance, Refs. 

While a big part of the theoretical investigations in this 
field deals with quantum systems in contact with bosonic 
environments, coherent transport induced by structured 
spin baths has for the first time been taken into con- 
sideration in Rcf. [3 . There, it was shown that, the 
transition probability from the lower to the upper level 
of a two-level system can be increased even for finite tem- 
peratures when modeling the environment as a series of 
spins 1/2 arranged in a star configuration. The authors 
also proved that the use of two separated baths, corre- 
lated with each other, can help to improve the transport 
compared to the use of a single environment. 

At the same time, it is known that decoherence is sen- 
sitive to quantum critical changes in the bath and that 
the correlation length of the environment can influence 
the rate of loss of quantumness in a system IB"21|. A 
possible justification for this behavior can be found in 
the monogamous nature of entanglement 22 1. 

Starting from the results of Rcf. 



step forward by considering the case of a spin bath that 
can experience a quantum phase transition and studying 
how a critical change in the environment will influence 
the noise-induced properties of a two-level system (TLS) 
coupled to the bath itself. The goal of our work is to find 
a direct link between the phase transition of the bath and 
the internal transition probability of the TLS. 



We approach the problem by first considering the case 
of a TLS coupled to a homogeneous XX chain in the pres- 
ence of a transverse field. If the system is coupled to the 
total magnetization of the bath, the problem is exactly 
solvable. We show that the transition probability inside 
the TLS is maximally sensitive to a change in the value 
of the applied field near its critical point and justify the 
result by noting that the bath renormalizes the system 
Hamiltonian parameters. The same qualitative behav- 
ior can be observed in the case of an inhomogeneous XY 
chain, where the analytical solution is no longer available, 
and where the simple mechanism of rcnormalization can- 
not be applied any more. 

Let us consider a TLS described by the Hamiltonian 



He = —a. 



(1) 



Because of the detuning, A, between the states | f) 
and | ],), a state initially prepared in | J,) has a maxi- 
mum probability of transition to | f) of P tr = T 2 /(T 2 + 
A 2 /4) < 1, with perfect transition for A = 0. As dis- 
cussed by Sinayskiy et al. in Rcf. jl5| . under certain 
circumstances, the presence of a single spin bath or two 
different, correlated baths, can help increase the maxi- 
mum transition probability for finite values of A. 

Here our aim is to analyze the effect an internally struc- 
tured environment can have on a decoherence-assisted 
transition process. We will consider a spin environment 
and study its influence as a function of an internal param- 
eter which drives the system through a quantum critical 
point. For this, we model the bath as an XY chain in the 
presence of a transverse field, which can be described as 
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where the total magnetization along z is M z 
and where the boundary conditions are periodic cr^ +1 = 
<7j*. From now on we will take J = 1 and use it as a 
scale of energy. In the thermodynamic limit the Hamil- 
tonian Hxy possesses a critical point which can be ob- 
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served when considering the average magnetization per 
spin, m z = M z /N H 



lim m z 
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(h — cos 0) d<^, (3) 



where /? 
[(ft- 



the inverse temperature 
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and A(/i) = 
For finite temperatures, the 



cos(/>)^ + 7^ sin" 
derivative of m z with respect to h shows a maximum for 
h < h Cl which turns into a divergence at h c = 1 for 
/3 — ^ oo. 

In the following, we will first pay special attention to 
the case of an isotropic chain, characterized by 7 = 0, and 
indicate the Hamiltonian operator as Hxx ■ For the sake 
of simplicity and without any loss of generality, we will 
only consider non-negative fields (h > 0) and note that 
M z is a conserved quantity since [Hxx, J2i a i) = 0- For 
h > h c = 1, the ground state is completely ordered along 
the direction z and no spin-spin correlations are present. 
Lowering h below the critical value h c , the ground state 
starts acquiring magnons (flipped spins) and for h = 
the number of magnons is exactly equal to N/2. In a pic- 
torial representation, exploiting the Jordan- Wigner map- 
ping between spins and fermions, the transverse field h 
plays the role of a chemical potential and determines the 
value of the Fermi level: all the levels with energy less 
than h will belong to the Fermi sea and will be filled, 
while all the levels with energy greater than h will be 
left empty. This mechanism implies that, for finite-size 
systems, moving the field from zero to one, N/2 tran- 
sition points between different symmetry sectors of the 
Hamiltonian are crossed. In the thermodynamic limit, 
this leads to an infinite-order quantum phase transition, 
named Berezinsky-Kosterlitz-Thouless, which takes place 
without spontaneous symmetry breaking. 

We assume that the interaction between the bath and 
the system is given by a coupling between the upper level 
t) and the total magnetization of the chain [15| 



flj = -r| t)(t I ® m. 



Since m z is a constant of motion, one can see that 
[Hxx, Hi] = and therefore the degrees of freedom of 
the environment can be eliminated. As a consequence, 
an exact solution for the dynamics of the system can be 
obtained. Notice that we are not dealing with a purely 
dephasing interaction, which would be characterized by 
[Hs,Hi] = 0, but that our model is suitable to describe 
a fully dissipative dynamical evolution. 

Assuming that at t = the system and the bath are 
in a product state [p{0) = ps(0) ® ps(0)], the density 
matrix of the system will evolve according to 



Ps(t) 



= Tr B [e-^ Hs+Hl)t p s (0) ® p B (0)e l{Hs+Hl)t ] 



(5) 



where we have now eliminated Hxx by exploiting 
[Hxx, HA = and the invariance of the trace under 



cyclic permutations. Even if formally eliminated from 
the dynamics, Hxx is present through the initial state of 
the bath ps(0), which will be assumed to be in a thermal 
distribution ps(0) = e~ l3Hxx /Z at the inverse tempera- 
ture /3, where the partition function is Z = Tr[e~^ Hxx ]. 

Both Hxx and Hi can be written in a diagonal form by 
using the aforementioned Jordan- Wigner transformation, 
which maps spins into spinless fermions [24j |. Labelling 
the eigenstates of the bath |A) and the corresponding 
cingcvalues A, the matrix elements (i\ps\j) 
evolve in time according to 



p l g will 
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- A, Vs(0)e^ A,t b-), (6) 



with Hg> =H s -T\ f)(t |(AK|A) = H s + T(2n x /N - 
1)| t)(t |j where n\ is the fermionic occupation number 
of |A>. 

Let us start analyzing the zero-temperature scenario, 
where the bath is in its ground state |G), which is char- 
acterized by its fermionic occupation number uq, which 
in turn depends on h. In this case, the sum over all 
the Hamiltonian eigenstates in Eq. ^ reduces to only 



one term, p l g{t) 
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This allows to 



by P tr = (i \e~ iH s ;t ps(0)e^ '*| t ) 

renormalisc the diagonal elements of the system Hamilto- 
nian and the effective gap between | t) and | l) becomes 
A = A - T(l - 2n G /N). For h = 0, where n G = N/2, 
the presence of the bath does not affect P tr , but for in- 
creasing field h the fermionic occupation number starts 
decreasing until we have tiq = at h = h c = 1. The 
value of n G /N in the thermodynamic limit for 7 — > 
and (3 — > 00 can be obtained from Eq. ((3]). If T > A it 
is therefore possible to achieve perfect energy transfer by 
choosing an optimal h between zero and one such that 
A = 0. For the case of T < A, which is the experimen- 
tally more realistic scenario, a maximum enhancement 
will be reached for h > h Cl which corresponds to the 
zero-temperature case described in Ref. [15| . 

For systems at finite temperatures we need to consider 
all terms of the sum in Eq. ©, which, for any finite N, 
can be evaluated analytically. For a finite inverse tem- 
perature (3 — 40, we show Pt r in Fig. [T]as a function of 
h for chains of 8, 12 and 16 spins and compare it to its 
value in the thermodynamic limit. The later is calculated 
by introducing a mean-field approximation, obtained by 
replacing Hi with H™ { = — T\ t)(t \ ( m z), where the aver- 
age magnetization per spin is given by Eq. ([3]) in the limit 
of vanishing 7. By comparing the different curves, it is 
evident that the curves for finite-sized chains are rapidly 
converging toward the mean-field approximation, which 
can then be used to estimate Pt r in the limit of long 
chains. The finite-size solutions show clear signatures of 
the N/2 transition points in the ground state between 
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FIG. 1: (Color online) P tr versus h for N = 8 sites [red 
(dashed) line], N = 12 sites [blue (dotted) line], N — 16 sites 
[green (dot-dashed) line], and for N = oo sites [black (solid) 
line], in the isotropic case (7 = 0). The system parameters 
are A = 2T = 0.2 and T = 5 x 10~ 2 , while the inverse bath 
temperature is j3 = 40. In the case of 8, 12, and 16 sites, Ptr 
has been calculated by using the analytical expression given 
in Eq. (JSJl , while the thermodynamic limit has been calculated 
within a mean- field approximation (see text for details). In 
the absence of the bath, one would have P tr — 0.5. Inset: 
the derivative of Ptr with respect to h in the thermodynamic 
limit. For finite-size chains, the derivatives would exhibit N/2 
peaks, the last of them falling around h c , in correspondence 
with the degeneracy points of the Hamiltonian. 



the different symmetry sectors of the Hamiltonian and 
exhibit horizontal plateaus between them. Taking the 
thermodynamic limit, where these transition points be- 
come infinitely dense, we are left with a critical change 
for Pt r at h ps h c . The existence of a critical point in Pt r 
is shown in the inset of Fig. [IJ where its derivative with 
respect to the external field is plotted. It shows a clear 
peak around h c , which is due to the behavior of m z . 

In the zero-temperature case, the behavior of P tr 
around h c is a direct consequence of the commutativ- 
ity between the bath and the interaction Hamiltonians. 
Indeed, the only effect of the bath is to replace the ini- 
tial, unperturbed, system Hamiltonian Hs with the effec- 

(G) 

tive two- level Hamiltonian H s , which has acquired the 
critical properties of the transverse magnetization. How- 
ever, in the finite-temperature case, this approach does 
not hold any more. Nevertheless, the results of Fig. Q] 
show that, for moderate temperatures, which can be ex- 
pected in realistic scenarios, a connection between P tr 
and the critical change in m z still holds (further decreas- 
ing /3 would lead the bath out of its quantum domain and 
wash out the effect on Pt r )- We will discuss this point in 
more detail below. 

Releasing the isotropy assumption by choosing 
7 7^ in Eq. ([2]), an analytical solution for 
Pg(t) analogous to that in Eq. ^ cannot be ob- 
tained. Instead, the general expression P tr — 
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FIG. 2: (Color online) Ptr versus h for N = 6 sites [blue (dot- 
ted) line], TV = 10 sites [red (dashed) line], and for N — 00 
sites [black (solid) line], in the presence of an anisotropic bath 
(7 = 0.8). The system parameters are A = 2T = 0.2 and 
F = 5 * 10~ 2 , while the inverse bath temperature is /3 = 40. 
Here, the finite-size lines have been calculated by explicitly 
solving the Hamiltonian evolution, while the mean-field ap- 
proximation has been used in the thermodynamic limit. Inset: 
Pt r as a function of 7 for strong fields (h — 2) and for /3 — 40. 
The transition probability is a monotonically decreasing func- 
tion of the anisotropy. 



max t \(i |Tr B [e-' iiftot Vs(0) <8> p B (0)e iHtatt }\ t}| 2 , where 
Htot = Hs + Hxy + Hi, must be evaluated. In the 
case of short chains, this can be done exactly by diago- 
nalizing -fftot and explicitly performing the trace over the 
environment. In Fig. [2] we show P tr for chains of 6 and 
10 spins at finite temperature and compare the results to 
the mean-field approximation, performed, as before, by 
replacing m z with its average value in Hi. One can see 
that the mean field solution is a good approximation to 
the exact, finite-size solution even for short chains, and 
therefore becomes a suitable candidate to describe the 
thermodynamic limit. Furthermore, this behavior shows 
that the influence of a quantum phase transition in the 
bath on the transition dynamics between the two internal 
system states is not limited to a special choice of Hamil- 
tonian or to possible special commutation properties, but 
lies in nature of the interaction. 

Let us finally discuss the interplay between tempera- 
ture effects and anisotropy for both small (h = 0.5) and 
large fields (h = 2). In Fig. [3] we show P tr as a function 
of j3 for the isotropic case 7 = and for the Ising case 
7 = 1. Two remarkable results emerge: i) the bath can 
have a positive effect on P tr for a large range of tem- 
peratures and ii) 7 is in fact important only in the very 
low temperature regime. The latter result not only holds 
for strong fields, where the x — y component of the in- 
teraction is dominated by the external field, but also in 
the symmetry-broken region, where the in-plane and the 
transverse terms have comparable strength. 

In conclusion, we have considered the problem of a 
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FIG. 3: (Color online) Pt r versus 1//3, in the thermodynamic 
limit, for 7 = and h — 2 [blue (dotted) line], 7=1 and 
h — 2 [green (dotdashed) line], 7 = and h = 0.5 [black 
(solid) line], and 7 = 1 and h = 0.5 [red (dashed) line]. 



TLS immersed in a spin chain that plays the role of the 
environment. We have shown that the transition proba- 
bility from the lower level to the upper level of the TLS 
is sensitive to the critical properties of the bath. De- 
scribing the bath as an isotropic XX chain and the in- 
teraction with the system through its transverse magne- 
tization, the problem becomes exactly solvable, since the 
bath Hamiltonian commutes with the interaction term. 
In the zero-temperature case, the bath renormalizes the 
system Hamiltonian making it possible to improve the 
transition probability. Even if the mechanism is more 
complicated, the same qualitative enhancement can be 
observed in the finite temperature regime. Finally, the 
decoherence-assisted process identified in our work has a 
more general range of validity, since we have shown that 
the same qualitative effect can be observed by using an 
anisotropic chain, where the commutativity property be- 
tween bath and interaction Hamiltonian does not apply 
any more. 

We believe that our results can elucidate the inter- 
play between quantum coherence and correlations and 
can help understand the basic mechanisms underlying 
energy transport in biologic systems. 
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